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ABSTRACT 

Pathogenic Escherichia coii strains carrying the 
afa-8 gene cluster are frequently associated with 
extra-intestinal infections in humans and animals. 
The afa-8 A to E genes determine the formation of 
an afimbrial adhesive sheath consisting of the AfaD- 
VIII invasin and the AfaE-VIII adhesin at the bacterial 
cell surface. This structure is thought to be required 
for host colonization. We characterized a new gene 
encoding the small RNA AfaR, which is transcribed 
in cis from the complementary strand of the 3' un- 
translated region of the afaD messenger RNA, 
within the afaD-afaE intercistronic region. AfaR is 
a frans-acting Hfq-dependent antisense small RNA 
that binds the 5' untranslated region of the afaD 
messenger RNA, initiating several ribonuclease 
E-dependent cleavages, thereby downregulating 
production of the AfaD-VIII invasin. AfaR transcrip- 
tion is dependent on a E , a member of the stress 
response family of extracytoplasmic alternative 
sigma factors. We found that the AfaR-dependent 
regulatory pathway was controlled by temperature, 
allowing the production of the AfaD-VIII invasin at 
temperatures above 37 C. Our findings suggest 
that the entry of afa-8-positive pathogenic E. coii 
strains into epithelial cells is tightly regulated by 
the AfaR small RNA. 

INTRODUCTION 

Escherichia coii has developed a commensal lifestyle in the 
lower part of the intestine of humans and other verte- 
brates, but some strains have the potential to cause a 
wide spectrum of intestinal and extra-intestinal diseases. 
Extra-intestinal pathogenic E. coii (ExPEC) strains cause 



bacteraemia, pyelonephritis, cystitis, prostatitis and 
neonatal meningitis in humans, and various infections in 
animals. The virulence of ExPEC strains is largely 
associated with the presence of virulence factors, including 
adhesins, toxins, siderophores, capsules, invasins and 
factors contributing to serum resistance (1,2). During the 
infection process, ExPEC strains must rapidly adapt the 
expression of their genes in response to both environmen- 
tal and host signals. A plethora of sensory systems for 
activating or repressing the expression of virulence genes 
has thus coevolved with virulence factors. 

In regulatory cascades, sensory systems and effectors 
may be connected by proteins. However, small non- 
coding RNAs (sRNAs) have recently emerged as a 
major class of gene regulators of adaptive responses, 
relaying information from physicochemical stress 
'sensing' components to response genes at the post- 
transcriptional level, without the need for translation (3). 
A significant proportion of the sRNAs characterized to 
date interact with dedicated mRNA targets through an 
antisense-based mechanism, affecting their stability and/ 
or translation (4). Each sRNA is thought to regulate the 
expression of more than one target, depending on its 
temporal pattern of expression, thus rendering the 
system even more complex. 

The genome-wide identification of new sRNAs in both 
gram-negative and gram-positive pathogenic bacteria, by 
biochemical and in silico approaches and, more recently, 
by high-throughput sequencing analyses, has highlighted 
the diversity and role of sRNAs (3,5-7). In particular, 
virulence-associated sRNA genes are abundant in the 
core genome and pathogenicity islands of the human 
pathogens Staphylococcus aureus (8-11), Streptococcus 
agalactiae (6), Streptococcus pyogenes (12), Salmonella 
typhimurium (13,14), Listeria monocytogenes (15), 
Clostridrium perfringens (16) and Vibrio cholerae (17). 
However, despite several investigations reporting the 
identification of new sRNAs in E. coii, only a few 
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virulence-associated sRNAs from pathogenic strains have 
been characterized to date. Studies of the role of the Hfq 
protein, which facilitates gene regulation by sRNAs, 
showed virulence to be strongly impaired in /^-deficient 
E. coli strains. These strains displayed defects in urinary 
tract colonization, motility and biofilm production [for 
uropathogenic E. coli UTI89 (18)], together with an im- 
pairment of epithelial cell invasion in the case of adherent- 
invasive E. coli LF82 (19). These data suggest that several 
sRNAs are involved in the control of virulence factor ex- 
pression by pathogenic E. coli/Shigella species. However, 
these molecules have been little studied, other than in 
quorum sensing (20), icsA -dependent invasion of the in- 
testine (21) and type 1-mediated adhesion (6). 

During host colonization, ExPEC strains produce 
several adhesins, including afmibrial adhesins of the Afa 
family (1,22). All afa gene clusters are organized in a 
similar manner, with six genes {afa A to F) (1) determining 
the formation at the cell surface of an adhesive sheath 
consisting of the AfaD invasin and the AfaE adhesin. 
These proteins are encoded by the afaD and afaE genes, 
and are responsible for bacterial internalization and 
binding to the host cell, respectively (23-26). Both these 
structural components are assembled via the chaperone- 
usher pathway, which is mediated by the proteins encoded 
by the afaB and afaC genes (27,28). Finally, the expression 
of afa gene clusters may be controlled by the afaF and 
afaA genes, which are homologues of the papI/papB tran- 
scriptional regulator genes (27,29,30). Various subtypes of 
the afa gene cluster (afa-1 to afa-S, dra, daa) have been 
described, on the basis of comparisons of the sequences of 
the afa gene cluster, focusing on afaE genes in particular 
(22). The afa-8 gene cluster is common in pathogenic afa- 
carrying strains isolated from humans and animals 
(30,31). ExPEC strains carrying afa-8 gene clusters are 
frequently detected in humans suffering from pyeloneph- 
ritis and septicaemia, but these clusters are absent from 
diarrhoea-associated strains (30). A recent screening for 
ExPEC-specific sRNA genes showed that some such 
genes were present in adhesin gene clusters, including the 
type 1 fimbria operon (6). In particular, a chromosome- 
borne sRNA candidate gene, SQ109, located in the afa-8 
gene cluster of the E. coli AL862 strain, has been identified 
as a putative riboregulator (6). 

We show here that SQ109 sRNA downregulates AfaD- 
VIII invasin production by acting as an Hfq-dependent 
antisense sRNA, inducing RNase E-dependent afaD 
mRNA decay. The expression of this sRNA gene is de- 
pendent on the sigma factor ct e and temperature, suggest- 
ing that AfaD-VIII production is controlled by the 
environment of the bacterial cell. Based on its location 
within the afa-8 gene cluster, we renamed this sRNA 
gene AfaR. 

MATERIALS AND METHODS 

Strains, plasmids and growth conditions 

The E. coli strains and plasmids used in this study are 
listed in Table 1. All strains were grown in Luria Bertani 
(LB) broth, at 37°C, with continuous shaking (140 rpm), 



and were harvested at the indicated OD 6 oo- We used the 
following antibiotics for plasmid selection: 100 ug/ml car- 
benicillin (Euromedex), 50 ug/ml kanamycin (Sigma), 
25 ug/ml chloramphenicol (Sigma), 12.5 ug/ml tetracycline 
(Sigma) and 100 ug/ml apramycin (Sigma). 

General DNA techniques 

Oligonucleotides were designed, and in silico analysis was 
carried out with the afa-8 sequence of E. coli 239KH89 
(accession number AF072900). The sequences of the oligo- 
nucleotides (Eurogentec) used in this study are listed in 
Supplementary Table SI. All polymerase chain reaction 
(PCR) products were amplified from the chromosomal 
DNA of E. coli strain 239KH89 (99.7% nucleotide 
sequence identity to strain AL862). Enzymatic reactions 
were performed in accordance with the manufacturer's 
instructions. 

Plasmids 

We generated the expression plasmid pILL1320 by amp- 
lifying the full afa-8 gene cluster by PCR with the primers 
clon.afa8.5 and clon.afa8.3 and Taq DNA polymerase 
(Qbiogen). The PCR product was digested with Hindlll 
(Roche), purified by agarose gel electrophoresis, and 
ligated into the Hindlll site of pBR322. 

For expression of the AfaR sRNA, a fragment contain- 
ing the afaR gene without its promoter was amplified by 
PCR with the Cl.afaR.EcoRI and Cl.afaR.Xbal oligo- 
nucleotides. The PCR product was purified by agarose 
gel electrophoresis and inserted into pCRII-TOPO 
(Invitrogen). The EcoRI/Xbal fragment of the pCRII- 
afaR plasmid containing the afaR gene was purified by 
agarose gel electrophoresis, ligated between the EcoRI 
and Xbal sites of pZE2R-g//? with T4 DNA ligase 
(Fermentas) and transferred into the TOP 10 strain by 
electroporation. 

We constructed the pZE2R-a/a^* plasmid, encoding 
the AfaR* mutated sRNA allele, as described for 
pZElR-afaR. We first introduced two mutations by amp- 
lifying the 5' and 3' regions of the afaR gene with the 
Cl.afaR.EcoRI/mutafaR5 and mutafaR3/Cl.afaR.XbaI 
oligonucleotides. The two resulting DNA fragments were 
assembled by PCR with the Cl.afaR.EcoRI and 
Cl.afaR.Xbal oligonucleotides. 

The pXGafaD::gfp fusion plasmid was constructed by 
amplifying the sequence of the intergenic region including 
parts of the 3' end of the afaC gene and the 5' end of the 
afaD gene with the afaC.Nsil and afaD.Nhel primers. The 
PCR product was inserted into pXG-30 digested with Nsil 
and Nhel, as previously described (34). The 
pXGafaD*::gfp fusion plasmid was constructed by 
introducing mutations into the afaD gene by the same 
strategy used for afaR*. 

The pXGlacZ::qfaD fusion plasmid was constructed in 
a similar manner. Briefly, pXG30 was digested with Nsil 
(Fermentas) and Xbal (Roche), and the 3' end of the afaD 
gene was amplified by PCR with the afaD. Nsil and 
afaD. Xbal oligonucleotides, treated with exonuclease I 
(Fermentas) and digested with Xbal (Roche). The 
digested plasmid and the PCR product were purified by 
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Table 1. Strains and plasmids used in this study 



Name 



Description/relevant characteristics 



Phenotype 



Source/reference 



E. coli strain 
239KH89 
ALIO 
AL213 
183 

AL511 

AL51 1 AafaR::KmFRT 

AL862 

N3433 

N3431 

BW25113 

JW4 130-1 
DH5ot 

TOP 10 



TOPl0Ahfq::KmFRT 
TOP\QAhfq::FRT 



MC1061 



CAG22216 
Plasmids 
pKD4 

pKOBEG-Apra 

pCP20 

pCP22 

pBR322 
pILL1320 

pILL1322 
pILL1323 



pILL1324 

pILL1325 

pZE2R-g#> 

pZE2R-null 

pZE2R-afaR 

pZE2R-afaR* 

pQF50 

pQFafaA::IacZ 

pQFafaB::lacZ 

pQFafaC.-.iacZ 

pQFafaD::lacZ 

pQFafaE::lacZ 

pQF afaEr.lacZ 

pQFafaR:.iacZ 

pXG-0 

pXG-30 

pXGafaD::gfp 

pXGlacZ ::afaD 



pXGafaD*::gfp 

pACYC184 
pTX381 



Bovine isolate, afa-8 + 
Human clinical isolate, afa-8 + 
Human clinical isolate, afa-8 
Human commensal isolate, afa-8 + 
Human clinical isolate, afa-8 + 
af a R-de&cient AL511 strain 
Human clinical isolate, afa-8 
lacZ43(Fs), LAM-, relAl, spoTl, thi-1 
N3433 me-3071(ts) 

K-12 A(araD-araB)567, lacZ4787A (::rrnB-3 j ,LAM-, rph-1, 

A(rhaD-rhaB), hsdR514 
BW25113 Ahfq-722::km 

F endAl hsdR17 supE44 thi-1 recAl gyrA relAl 

AflacIZYA-argF)U169 deoR {<&&0/acZAM]5) ().pyr + ) 

mcrA A( mrr-hsdRMS-mcrBC) Cp80/flfZAM15 AlacXIA 
deoR recAl araD139 A (ara-leu)l 697 galU galK rpsL 
endAl nupG 

///^-deficient TOP 10 strain 

///^-deficient strain with the FRT-flanked kanamycin 
cassette removed by flippase-mediated excision with 
pCP20 

K-12 F- araD139 A(ara-leu)7697 A(codB-lacI) galK16 
galE15 mcrAO relAl rpsL150 spoTl mcrB9999 hsdR2 
lacX74 

MC1061 [fykrpoH P3::lacZ\ ArpoEr.Cm 

Source of the KmFRT resistance cassette 

Derivative of pKOBEG 

FLP flippase expression plasmid 

pCP20 plasmid with insertion of the ble gene into Xhol site 

conferring resistance to zeocin 
Cloning plasmid 

afa-8 gene cluster from 239KH89 inserted into Hindlll site 

of pBR322 
pILL1320 ARx-afaR::FRT 

pILL1320 AafaR::FRT, complete deletion of the region fol- 
lowing the rho-independent terminator of afaD gene with 
partial removal of afaR gene sequence 

pILL1320 AafaD::FRT 

pILL1322 ARx-afaA::FRT 

gfp gene under control of the PA constitutive promoter 
plasmid encoding a ~60nt nonsense RNA 
afaR gene under control of the PA constitutive promoter 
p7.F2R-afaR derivative plasmid carrying the GT point 

mutated allele of AfaR 
LacZ fusion plasmid 

5' UTR of afaA inserted into Xhol/Xbal sites of pQF50 
5' UTR of afaB inserted into Xhol/Xbal sites of pQF50 
5' UTR of afaC inserted into Xhol/Xbal sites of pQF50 
5' UTR of afaD inserted into Xhol/Xbal sites of pQF50 
5' UTR of afaF inserted into Xhol/Xbal sites of pQF50 
5' UTR of afaE inserted into Xhol/Xbal sites of pQF50 
5' region of afaR inserted into Xhol/Xbal sites of pQF50 
Translational fusion control plasmid, gfp~ 
Translational fusion plasmid, gfp + 

Insertion of the 5' UTR part of afaD into Nhel/BfrBI sites 
of pXG30, giving a afaDr.gfp translational fusion 

Insertion of the 3' UTR part of afaD into Nhel/Xbal sites 
of pXG30, giving a M2-tagged lacZr.afaD translational 
fusion 

pXGafaDr.gfp derivative plasmid carrying the AC point 

mutated allele of AfaD 
Cloning plasmid 

hfq gene inserted into pACYC184 



Sm R 



Km" 
Km s 



Cm R 
Km R 

Apra R , Cm R 
Cb R 

Cb R , Zeo R 
Cb R , Tet R 



Cb R , Tet s 
Cb R , Tet s 



Cb R , Tef 

Cb R , Tet s 

Km R 

Km R 

Km R 

Km R 



Cb K 
Cb R 
Cb R 
Cb R 
Cb R 
Cb R 
Cb R 
Cb R 
Cm R 
Cm R 
Cm" 



Cm" 



Cm" 



Tet R , Cm R 
Tet R , Cm R 



29 
31 
31 
32 
31 

This study 
29 

CGSC# 6976 
CGSC# 6975 
33 

33 

Lab. collection 
34 

34 

This study 
Lab. collection 

3? 

36 
37 
38 

This study 

Lab. Collection 
This study 

This study 
This study 



This study 
This study 
39 
6 

This study 
This study 

40 

This study 
This study 
This study 
This study 
This study 
This study 
This study 
34 
34 

This study 
This study 



This study 

Lab. Collection 
41 



CGSC, Coli Genetic Stock Center. 

Phenotypes: Apra, apramycin; Cb, carbenicillin; Cm, chloramphenicol; Km, kanamycin; Sm, streptomycin; Tet, tetracycline; Zeo, zeocin. S, sensitive; 
R, resistant. 
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agarose gel electrophoresis, ligated with T4 DNA ligase 
(Roche), and introduced into the TOP 10 strain by 
electroporation. 

All the promoter fusion plasmids constructed were de- 
rivatives of pQF50. The 300 nucleotides upstream from 
the initiation codon of the afaA to E coding sequences 
(CDSs) and the promoter region of afaR were amplified 
by PCR with the primers listed in Supplementary Table 
SI. PCR products and pQF50 were digested with Xhol 
and Xbal, purified by agarose gel electrophoresis, ligated 
with T4 DNA ligase (NEB) and used to transform E. coli 
strain DH5a by electroporation. 

All plasmid constructs were sequenced (GATC Biotech) 
to check that there were no unwanted mutations. 

Promoter fusion plasmid assay 

We analysed B-galactosidase activity on LB agar plates 
containing 40ug/ml X-Gal and lOOug/ml carbenicilin by 
scanning bacterial colonies with a PharosFX (Biorad) at 
an excitation wavelength of 480 nm. The E. coli DH5oc 
strain bearing pQF50 was used as a negative control. 

Inactivation of genes by allelic exchange 

Mutant strains were constructed by the allelic exchange 
recombination method, with the thermosensitive plasmid 
pKOBEG-Apra. Briefly, long synthetic primers (~60nts) 
overlapping the 5' and 3' untranslated regions of target 
genes were used to amplify the kanamycin FRT-flanked 
cassette (KmFRT) from pKD4 by PCR. Flippase- 
mediated excision of the kanamycin resistance cassette 
was achieved with pCP22. 

RNA isolation and reverse transcription 

Total RNA was extracted by the hot-phenol protocol (6) 
and treated twice with DNase I (Roche). We checked the 
quality of total RNA samples on a Bioanalyzer with an 
RNA 6000 nanochip (Agilent). Genomic DNA contamin- 
ation was assessed by PCR amplification of the 5S ribo- 
somal gene with the 5S.Fw and 5S.RT primers, and was 
considered insignificant if no PCR product was observed 
on agarose gel electrophoresis. Chimeric DNAs were 
synthesized from total RNA with random primers for 
mRNA and specific primers for sRNA, as previously 
described (6). 

Quantitative RT-PCR 

For quantitative reverse transcription PCR (qRT-PCR) 
analyses of gene expression, we used Primer3 software to 
design specific primer pairs. Real-time PCR was per- 
formed with the Bio-Rad iQ system and iQ SYBR 
Green supermix (Bio-Rad), according to the manufac- 
turer's instructions. The thermocycling protocol was as 
follows: denaturation at 94°C for lOmin, 40 cycles of 
amplification and quantification (95°C for 15 s, 60°C for 
1 min with a single fluorescence measurement) and a 
melting curve program (55°C to 95°C, in increments of 
0.5°C). Each assay was performed at least in duplicate, 
with two independently prepared total RNA samples. 
The relative expression data obtained were analysed with 



the optimized amplification efficiency protocol derived 
from the AACt method, with the 5S gene used as a refer- 
ence gene (42). 

5' and circular RACE 

The 5' RACE experiments were carried out as previously 
described, but with the following modifications (43). 
Briefly, 10 ug of freshly prepared total RNA from the 
239KH89 strain was left untreated or treated with the 
TAP enzyme (Epicentre) and ligated with the Rl RNA 
adapter (Sigma Proligo). Chimeric DNA was synthesized 
with Superscript 111 reverse transcriptase (Invitrogen) and 
random primers (Invitrogen). It was then amplified by 
PCR with P10 and reverse primers and analysed by 
agarose gel electrophoresis. The largest TAP-treated 
PCR products were extracted from the gel, inserted into 
pCRII-TOPO (Invitrogen) and sequenced (GATC 
Biotech). 

Circular RACE was carried out on freshly prepared 
total RNA. The RNA was first denatured by heating at 
80°C for 2 min in lx buffer (10 mM Tris-HCl pH 8.1 at 
37°C, 7% deionized formamide, 0.2% SDS) and was then 
dephosphorylated with CIP enzyme (NEB). The RNA was 
purified by phenol/chloroform extraction and isopro- 
panol precipitation. The 5' ends of RNAs were 
re-phosphorylated with T4 polynucleotide kinase (NEB) 
in lx buffer (80 mM Tris-HCl pH 8.1 at 37°C, 8mM 
MgCl 2 , 4mM dithiothreitol, 1 mM spermidine), 1 mM 
ATP and 40 units of RNaseOUT (Invitrogen). The 
RNA sample was self-ligated, with or without T4 RNA 
ligase (Epicentre). Circular RNA products were reverse 
transcribed with random primers, and the products were 
amplified by PCR with rigBC.RT and rigCD.Fw primers 
and checked by agarose gel electrophoresis. The DNA of 
bands of interest (detected only after the use of T4 RNA 
ligase) was extracted from the agarose gel with the Qiagen 
gel extraction kit, inserted into pCRII-TOPO and 
sequenced. 

Measurement of RNA decay by rifampicin assay 

We measured the rate of decay of AfaR sRNA by growing 
strain carrying the pZE2R-afaR plasmid in LB broth, at 
37°C or 42°C, with continuous shaking (140 rpm) until the 
OD 60 o reached 0.4. Rifampicin was added to stop the ini- 
tiation of transcription (final concentration, 500ug/ml). 
Aliquots of cells were removed at several time points 
(0, 1, 2, 4, 8, 16, 30 and 60 min) after the addition of 
rifampicin and immediately frozen by mixing with an 
equal volume of cooled absolute ethanol. RNA prepar- 
ation, cDNA synthesis and quantitative RT-PCR were 
performed as described above. The half-life was calculated 
from the slope of a least-squares regression line of a semi- 
logarithmic plot of the percentage sRNA remaining as a 
function of time. Half-life was calculated from the time 
points at which decay rate was exponential. 

Western blotting 

Crude extracts of bacterial cell were separated by 
SDS-PAGE in 14% polyacrylamide gels, and the 
protein bands were electro-transferred onto PVDF 
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membranes (Millipore). The membranes were then 
hybridized with the a-GFP antibody (Roche), the oc- 
FLAG antibody (Sigma-Aldrich) and oc-GroEL antisera 
(Sigma-Aldrich) as previously described (34). The AfaD- 
VIII and AfaE-VIII proteins were detected by incubation 
with rabbit polyclonal anti-rAfaD-VIII or anti-rAfaE- 
VIII antibodies (25) at a dilution of 1:2000, followed by 
horseradish peroxidase (HRP)-conjugated anti-rabbit 
serum (Zymed). The membranes were then developed 
with the SuperSignal West pico chemiluminescent 
reagent (Pierce), and the signals were detected with a 
ChemiDoc XRS system (Biorad). 

RESULTS 

Transcriptional organization of the afa-8 gene cluster 

The sequence and genetic organization of the afa-8 gene 
cluster of E. coli 239KH89 strain were described a few 
years ago, but we still know little about its transcription 
(29). We first predicted the transcriptional units of this 
afa-8 gene cluster in silico by determining the location of 
the promoters and Rho-independent terminators. Putative 
a 70 promoters were identified upstream from the afaA, 



afaE and afaF genes (Figure 1A). These promoters were 
validated by fusing the 300-bp DNA fragment upstream 
from the AUG start codon of each afa CDS with the lacZ 
reporter gene of the pQF50 plasmid. (3-galactosidase 
activity was observed only for strains carrying 
pQFqfaA:.iacZ, pQFafaE::lacZ and pQFafaF::lacZ 
(Figure IB), confirming our predictions. The identification 
of a Rho-independent terminator beginning 240 bp down- 
stream from the stop codon of the afaD gene suggests that 
the afaD mRNA has a long 3' UTR. Another Rho-inde- 
pendent terminator was found 25 bp downstream from the 
afaE stop codon. In a previous study, the ~230-nt SQ109 
sRNA gene was identified between the afaD and the afaE 
genes (6). This sRNA was renamed AfaR (afimbrial 
adhesin small RNA), as its location suggested a 
probable role in regulating the expression of the afa-8 
gene cluster. BLASTN analysis against the GenBank 
database (24 January 2012) revealed that the sequence 
and synteny of the afaR gene were conserved in all 
sequenced afa-8 gene clusters from E. coli isolates 
(239KH89, AL511 and AL862 strains) but that this gene 
was absent from all other afa family gene clusters. We 
mapped the transcription start site of afaR by 5' RACE 
to nucleotide position A5534 of the afa-8 gene cluster of 




Figure 1. Organization of the transcriptional units of the afa-8 gene cluster from E. coli 239KH89. (A) Genetic organization of the afa-8 gene cluster, 
with predicted promoters (arrows) and rho-independent terminators (circles). Position of primers (vertical dashed line) for PCR products [horizontal 
black lines 1 to 7 refer to the RT-PCR shown in (D) and (E)]. (B) p-Galactosidase plate activity assays. Each MC1061 strain carries a pQF50-derived 
plasmid with or without the region upstream from an afa gene fused to lacZ. (C) 5' RACE for identification of the AfaR transcription start site. 
Total RNA from E. coli strain 239KH89 was linked to a 5' RNA adaptor with [TAP (+) lane] or without [TAP (-) lane] removal of the 5' 
triphosphate by TAP treatment. After RT-PCR amplification with specific primers, PCR products were analysed by electrophoresis in 4% agarose 
gels. A single ~95 bp band strongly enhanced by TAP treatment was cloned and sequenced to identify the transcription start site of the afaR gene. 
(D, E) Cotranscription analysis of the afaA to E genes from the afa-8 cluster. (D) Analysis of the BW25113 + pILL1320 strain. (E) Analysis of the 
BW25113 + pILL1322 strain. The products of RT-PCR amplification were analysed by electrophoresis in a 1% agarose gel, which demonstrated 
probable linkage between afaA and afaB (lane 1), afaB and afaC (lane 2), afaC and afaD (lane 3), afaD and the 3' UTR of afaD (lane 4), the 3' UTR 
of afaD and afaE (lane 5), afaA and afaC (lane 6), afaA and afaD (lane7). The expected PCR products are indicated by boxes. 
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E. coli 239KH89 (Figure 1C). An analysis of the afaR 
promoter region showed strong sequence identity to the 
consensus sequence of the cj e (also named a 24 ) promoter. 
In particular, this region contains the 'AA' tract of 
the -35 element, the invariable C-residue of the -10 box, 
the 16-bp spacer, the 6-bp discriminator sequence and the 
conserved -1 T-residue (44). The AfaR promoter was 
validated by constructing the pQF afaR::lacZ fusion 
plasmid, which yielded P-galactosidase expression 
(Figure IB). The distance between the a E promoter and 
the Rho-independent terminator confirmed the previous 
estimate of 230 nt for the size of the AfaR sRNA (6). 

Owing to its antibiotic multiresistance phenotype, 
genetic modification of E. coli strain 239KH89 was not 
possible. We therefore carried out most of the experiments 
on the E. coli BW25113 laboratory strain carrying 
pILL1320, corresponding to the entire afa-8 gene cluster 
inserted into the pBR322 plasmid. Several RT-PCRs were 
performed to determine whether the afaABCDE genes 
were cotranscribed in the BW25113 + pILL1320 strain; 
for these reactions, we used primers targeting the 5' and 
3' ends of each gene and the 3' UTR of afaD (Figure 1A 
and D). The specificity of the reverse transcription reac- 
tions was demonstrated in control experiments by 
omitting the reverse transcriptase enzyme (data not 
shown). We checked that the transcription profile of the 
afaABCDE genes was similar in the E. coli 239KH89 and 
BW25113 + pILL1320 strains (Supplementary Figure SI). 
Our findings suggested that afaA and afaB, and afaB and 
afaC were cotranscribed, generating a single mRNA 
molecule. However, no linkage was detected between 
afaC and afaD, afaD and afaE or between afaA and 
afaD. Amplification of the afaD and afaD 3' UTR 
mRNAs but not of the afaD 3' UTR and afaE mRNAs 
confirmed the location of the Rho-independent terminator 
and suggested that afaE was independently transcribed. 
qRT-PCR analysis of the transcription of the afa-8 gene 
cluster of the BW25113 + pILL1320 strain revealed that 
relative transcript levels for the afaAB, afaC and afaD 
units were significantly different (Supplementary Figure 
SI). We interpreted these observations as indicating that 
a putative afaABCD mRNA could be cleaved by one or 
more endoribonucleases (RNases) at sites surrounding the 
afaC mRNA, resulting in the release of the afaD mRNA. 
The cotranscription of the afaABCD mRNA did not rule 
out the presence of a putative internal promoter. We 
therefore deleted the entire promoter region between 
afaE and afaA and the region downstream from the 
afaD Rho-independent terminator, including the afaR 
gene promoter from the pill 1320 plasmid, giving 
pILL1325. Analysis, by qRT-PCR, of the transcription 
of the afa-8 cluster genes from the BW25131 +pILL1325 
revealed no significant amplification of afaABCD genes 
with respect to the 5S gene, indicating that there was not 
such internal promoter in the gene cluster. In addition to, 
the afaE gene was amplified, confirming that it was an 
independent gene (data not shown). Taken together, 
these results suggest that the afa-8 gene cluster consists 
of four transcriptional units: afaABCD, afaR, afaE and 
afaF. 
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Figure 2. Predicted RNA secondary structures. (A) Nucleotide pairing 
of AfaR with the 5' UTR of the afaD mRNA from E. coli 239KH89. 
The putative RBS required for afaD translation and the ATG codon are 
highlighted in grey. Numbers denote the residues in the afaD mRNA 
with respect to the ATG codon and the mapped +1 site of AfaR. The 
positions of point mutations introduced into the afaD gene (TG 480 3 to 
AC, giving the afaD* allele) and AfaR (AC 110 to TG, giving the AfaR* 
allele) and expected to maintain base-pairing between the n/aD*/AfaR* 
duplex, are indicated by thin black arrows. The light grey arrow indicates 
the putative RNase E cleavage site at G 4791 (-32 relative to the ATG 
codon). The AU-rich sequence next to a stem-loop helix and a G residue 
positioned two nucleotides upstream from a cleavage site mimicking the 
RNase E recognition sequence are displayed in italics. The G residue is 
shown in italics and underlined. The stop codon of the afaC gene is 
underlined. (B) Predicted secondary structures of AfaR sRNA. The 
region of the afaD mRNA to which AfaR binds is underlined. 
Covariations supporting structure prediction are highlighted in grey. 

AfaR interacts with the afaD mRNA via an antisense 
mechanism 

Our working hypothesis is that AfaR is an antisense RNA 
that regulates AfaD production. We therefore used 
intaRNA software (45) to analyse in silico the likelihood 
of AfaR forming double-stranded antisense pairs with 
other afa-8 gene mRNAs. One significant putative antisense 
pairing was predicted between nucleotides U 92 to Ci 30 of 
AfaR and nucleotides G_ 37 to A +1 of the afaD mRNA 5' 
UTR, which contains the ribosome-binding site (RBS) 
(Figure 2A). A second predicted pairing involved the 3' 
UTR of the afaD mRNA (residues T 490 to G 6 95), which 
perfectly matched the last 206 nt of the AfaR sequence 
(data not shown). This putative pairing was predicted on 
the basis of the presence of afaD and the 3' end of the afaR 
gene in the same DNA locus, but in opposite orientations. 
Both pairing regions were conserved in E. coli strains 
carrying afa-8 gene clusters, as were AfaR secondary struc- 
tures (Figure 2B). We thus hypothesized that AfaR bound 
the 5' and/or 3' UTR of the afaD mRNA, thereby 
modulating its translation and/or stability. 
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For the validation of these predictions in vivo, we 
assessed the putative base-pairing interaction between 
AfaR and these two mRNA targets, using a translational 
control and target recognition system (34). The afaD 
5' region, spanning the last 66 bp of the afaC gene and 
the first 60 bp of the afaD gene, was inserted into pXG- 
30, giving pXGafaDrgfp. The expression of the afaDr.gfp 
translational fusion gene was under the control of the 
Pueto-i promoter. The afaR gene was placed under 
the control of the strong constitutive P k promoter in the 
pZE2R-afaR plasmid. Four E. coli TOP 10 strains 
harbouring a combination of the pXGafaDrgfp target 
plasmid or pXG-0 (no target control) and either pZE2R- 
afaR or pZE2R-null were constructed. We monitored the 
levels of the sRNA and fusion mRNAs by qRT-PCR and 
GFP production by Western blotting (Figure 3A). AfaR 
over-expression was associated with a decrease in the 
amount of the afaDr.gfp fusion mRNA by a factor of 2 
to 2.56 (Figure 3A). Interestingly, AfaR levels decreased 
with those of the afaD mRNA target, suggesting 
AfaR-guided ribonuclease degradation of the target. 
Western blot experiments with antibodies directed 
against GFP showed that AfaR overexpression was 
associated with a halving of AfaD::GFP protein levels, 
consistent with the results of the transcriptional analysis 
(Figure 3A). 

The AfaR/afaD RNA interaction was validated by 
introducing point mutations into the pairing regions of 
the AfaR gene and the afaD mRNA 5' UTR 
(Figure 2A). The TG_j 9 dinucleotide (UG in RNA) of 
afaD was changed to AC, yielding the afaD*rgfp allele, 
and the CA 10 9 dinucleotide in afaR was changed to GT, 
giving a compensatory afaR* allele. Western blot analyses 
showed that regulation of the afaDr.gfp fusion by AfaR* 
was strongly impaired. Similarly, the afaD*::gfp fusion 
was more resistant to regulation by wild-type AfaR 
allele expression from pZE2R-afaR. The regulation of 
afaD*::gfp fusion expression was restored by expressing 
the compensatory AfaR* allele from pZE2R-afaR* 
(Figure 3B). 

The 3' UTR of the afaD mRNA was the second 
putative target of AfaR determined in silico. We 
investigated whether the afaD 3' UTR played a critical 
role in regulation by inserting the sequence spanning the 
last 60 bp of the afaD gene and the following 278 bp 
sequence, including the afaD rho-independent terminator, 
into pXG30, to obtain pXGlacZrafaD. Expression of the 
lacZr.afaD fusion was monitored by qRT-PCR and 
Western blotting in E. coli TOP 10 harbouring a combin- 
ation of the pXGlacZ::qfaD target plasmid or pXG-0 
(control containing no target sequence) and either 
pZE2R-afaR or pZE2R-null (Figure 3C). A comparison 
of the relative levels of lacZr.afaD mRNA in strains 
bearing pZE2R-afaR or pZE2R-null showed that AfaR 
overexpression had no effect on fusion mRNA levels 
(Figure 3C). Western blots with antibodies directed 
against the M2-FLAG tag (LacZ) confirmed this observa- 
tion (Figure 3C). Thus, AfaR is an antisense sRNA that 
represses afaD expression in vivo by interacting with the 5' 
UTR of the afaD mRNA but not with its 3' UTR. 



AfaR controls afaD mRNA levels posttranscriptionally 

We constructed an afaR mutant strain with a view to con- 
firming the target of AfaR in vivo. E. coli strain 239KH89 
is resistant to several antibiotics, AL511 does not express 
afaE (Supplementary Figure SI) and AL862 carries two 
afa-8 gene clusters (29). We therefore mutated the 
sequence inserted into pILL1320. As the 3' sequences of 
the afaD and afaR genes overlapped, we had to delete the 
cr E promoter of afaR to abolish the expression of this 
gene, to generate pILL1322 (Figure 4A). We checked 
that the afaR sequence in the pILL1322 plasmid was not 
expressed by qRT-PCR (Figure 4B). A transcriptional 
linkage analysis of BW25113 + pILL1322 showed that 
afaC and afaD were cotranscribed. This suggested that 
AfaR might be involved in the release of afaD mRNA 
from a primary transcript (Figure IE), consistent with 
an interaction between the two RNAs based on an 
antisense mechanism. We then determined the relative 
levels of expression of the afaABCDE genes in 
BW25113 + pILL1320 and BW25113 + pILL1322 by 
qRT-PCR (Figure 4B). RNA was extracted from 
cultures at an OD 600 of 0.6, in LB medium, in which 
AfaR levels were maximal (data not shown). No signifi- 
cant difference was observed in the relative levels of ex- 
pression of afaABC and afaE, indicating that AfaR had 
little transcriptional or posttranscriptional effect on the 
expression of these genes. 

However, the inhibition of AfaR increased afaD expres- 
sion by a factor of 3.8 (Figure 4B), and Western blotting 
indicated an accumulation of the AfaD protein (~ 1 5 kDa) 
by a factor of 2.6 in the BW25113 + pILL1322 strain 
(Figure 4C). We validated these results by constructing 
the AL511 AafaRrKmFRT strain and analysing the ex- 
pression of the afa-8 gene cluster of this mutant by qRT- 
PCR, comparing the results obtained with those for the 
wild-type AL511 strain. The transcription levels of the 
afaA, B and C genes were similar in the two strains. By 
contrast, the afaD mRNA was overproduced in the afaR- 
lacking AL511 strain, as previously reported for 
BW25113 + pILL1322 (Figure 4D). Taken together, 
these results suggest that AfaD expression is controlled 
posttranscriptionally by AfaR. 

The pairing of the afaR and afaD mRNAs is dependent 
on Hfq 

In E. coli, most antisense sRNAs make use of the Hfq 
protein to capture their mRNA targets (46). We assessed 
the contribution of Hfq to the AfaR-mediated regulation 
of afaD mRNA in vivo by comparing the levels of afaD 
transcripts in the /^-depleted strain JW4130-1 and the 
isogenic BW25113 strain, both carrying pILL1320. The 
relative level of afaA mRNA in the hfq mutant was one 
tenth that in the wild type, as shown by qRT-PCR, 
indicating that Hfq affected the transcription of the afa- 
8 gene cluster (data not shown). When the pTX381 
plasmid (carrying the E. coli hfq gene) was introduced 
into the JW4130-1 +pILL1320 strain, afaA mRNA 
levels were restored to wild-type values (data not 
shown). We therefore normalized gene expression with 
respect to afaA in qRT-PCR analyses, to evaluate the 



5476 Nucleic Acids Research, 2013, Vol. 41, No. 10 



afaDr.gfp 



a-GEP o 
\ 

a-GroEL ^ 



- afaDrgfp 



B TOP10 

afaD::gfp afaD i: ::gfp 



11 

0.0 -I— 



:afaD 



11 



±t ri- m m 58 



Figure 3. AfaR-induced afaD mRNA decay occurred owing to the 
binding of AfaR to the 5' UTR of the afaD mRNA. (A) Analysis of 
binding to the 5' UTR of afaD mRNA. Analysis, by qRT-PCR and 
Western blotting, of afaDr.gfp and AfaR gene expression in E. coli 
strain TOP10 harbouring pZE2R-afaR or pZE2R-null together with 
pXG-0 (no target control) or pXGafaD::gfp target expression 
plasmids. The four isolates were cultured in LB medium at 37°C, to 
an OD 60 o of 1.0. The expression of the gfp fusion gene was normalized 
to 1.0 for the TOP10 vector carrying pZE2R-null and pXGafaD::gfp. 
AfaR expression was normalized to 1.0 for TOP10 carrying pZE2R- 
afaR and pXG-0. The GroEL protein was used as a loading control. 
(B) Western-blot analysis with antibodies directed against GFP, for 
TOP10 harbouring pTE2K-afaR or mutated pZE2K-afaR* together 
with either wild-type pXGafaD::gfp or mutant pXG afaD*:: gfp fusion 
plasmids. (C) Analysis of binding to the 3' UTR of afaD mRNA. qRT- 
PCR and Western blot analysis were performed as described in (A), but 
with strains carrying pXGlacZ::afaD instead of the pXQafaD::gfp 
fusion. 



transcription-independent effects of Hfq (Figure 5A). The 
deletion of hfq had a significant effect on the relative levels 
of afaC and afaD mRNAs and AfaR (Figure 5A), sug- 
gesting possible cooperation between AfaR and Hfq. We 
tested this hypothesis by carrying out the same experiment 
with the JW4130-1 + pILL1322 strain, in which expression 
levels were similar to those for JW4130-1 +pILL1320 
(Figure 5B). We characterized this effect in more detail 
by transforming the TOP 10 Ahfq strain harbouring 
either pZE2R-afaR or pZE2R-null with pXGafaD::gfp 
and monitoring the expression of the afaDr.gfp fusions 
and AfaR by qRT-PCR. The overtranscription of AfaR 
did not significantly decrease the amount of fusion 
mRNA. Furthermore, AfaD::GFP fusion protein levels 
were slightly higher in the presence of the AfaR sRNA, 
as shown by Western blots with a-GFP serum. The spot 
intensities for the fusion proteins in each lane were 
normalized with respect to the GroEL loading control 
for the lane concerned. The differences in normalized 
spot intensities were not statistically significant (Mest) 
(Figure 5C). We infer from these findings that the Hfq 
protein is required for the AfaR-mediated regulation of 
afaD expression. 



AfaR and RNase E-dependent cleavages of afaD mRNA 

Characterization of the transcriptional units of the afa-8 
gene cluster suggested a possible role for RNases in the 
posttranscriptional regulation of afaABCD mRNA 
(Figure 1C and D). We tested this hypothesis by analysing 
the dependence of afa-8 cluster gene expression on RNase 
E, an RNase frequently associated with sRNA mechanisms. 
We introduced pILL1320 into the N3431 strain carrying the 
thermosensitive allele me-3071 (encoding RNaseE ts ) and the 
isogenic strain (N3433) carrying the wild-type rne gene 
(encoding RNase E). Both strains were grown to an 
OD 60 o of 0.6 at 30°C, and half the culture was then 
shifted to 42° C for 30min to inactivate rne ts gene expres- 
sion. We analysed the relative levels of afaA and afaD 
mRNAs and of AfaR in the two strains by qRT-PCR. In 
the control RNase E + strain context (Figure 6C), afaA ex- 
pression levels were unaffecting, indicating that a shift in 
temperature to 42° C had no visible effect on transcription 
of the afaABCD operon. By contrast, this shift had a strong 
effect on AfaR levels. When the rne ts gene was inactivated, 
we observed a small, but marginally significant, increase by 
a factor of 2.4 in the relative level of the afaA transcript 
(Figure 6B). By contrast, the levels of the afaD and AfaR 
RNAs were 4.6 and 16.8 times higher, respectively 
(Figure 6B), indicating a possible role of RNase E in the 
posttranscriptional processing of afaD mRNA. A similar 
analysis was carried out with the RNase E ts strain 
carrying the pILL1322 plasmid. The levels of the afaA 
and afaD mRNAs did not change significantly with 
temperature, indicating that AfaR was required for 
RNase E-dependent cleavage of the afaABCD mRNA 
(Figure 6D). 

We assessed the possible AfaR-dependent cleavage of 
the afaABCD mRNA by carrying out differential circular 
RACE experiments between strains BW25113 carrying 
pILL1320 or pILL1323, a derivative of pILL1320 contain- 
ing only the complete afaABCD operon with its natural 
rho-independent terminator (Figure 6A). The use of this 
construct overcomes problems due to the influence of 
downstream sequences on the non-specific expression of 
AfaR. RNase cleavage occurred around afaC. We there- 
fore grew both isolates to an OD 600 of 0.6, extracted total 
RNA and analysed both AfaR levels and the ligation-de- 
pendent circularization of the afaC mRNA (See Methods 
section, Figure 6E). We sequenced 20 recombinant 
plasmids containing the ligated 5' and 3' ends of the 
afaC mRNA (Supplementary Figure S2). We filtered out 
sequences that were not identified at least three times, to 
distinguish between AfaR-dependent RNase activity and 
non-specific cleavage. Three cleavage sites were identified 
at positions U_278o, G_32 and A 66 , within the 5' end of the 
afaC mRNA, the 5' end of the afaD mRNA that paired 
with AfaR and within the coding region of the afaD 
mRNA, respectively (Figure 4A and 6E). Relative levels 
of afaB and afaC mRNA in BW25113 + pILL1322 and 
BW25113 + pILL1320 were similar, indicating that 
cleavage in the 5' UTR region of the afaC mRNA was 
independent of AfaR (Figure 4B). For confirmation that 
afaD mRNA cleavage was dependent on AfaR, we 
performed a differential 5' RACE experiment on 
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Figure 4. Posttranscriptional control of afaD mRNA levels by AfaR, in vivo. (A) Overview of the afaDjafaR locus from the afa-8 gene cluster of 
E. coli strain 239KH89. The promoter boxes (-10 and -35) and transcription start site (+1) of the afaR gene are underlined. Rho-independent 
terminators are indicated by an arrow. 'ATG' in bold and asterisks indicate the initiation and termination codons, respectively, of the afaC and afaD 
CDSs. The grey box indicates the binding region of the AfaR sRNA and afaD mRNA. Vertical grey arrows indicate the position of ribonuclease 
cleavages. The pILL1322 plasmid has the same sequence as pILL1320, except that the a E promoter sequence is replaced by the FRT site, as 
indicated. (B) Assessment of the relative expression of the afaABCDER genes in strain BW25113 carrying pILL1322 (APr-afaR) versus pILL1320 
(wild type), determined by qRT-PCR. Bacteria were grown to an OD 60 o of 0.6. (C) Western blot analysis (with an anti-rAfaD-VIII antiserum) of 
AfaD production in a total protein extract obtained at an OD 600 of 0.6 from BW25115 carrying pBR322 (mock), pILL1320 (wild type), pILL1322 
(AVr-afaR) or pILL1324 (AafaD). The GroEL protein was used as a loading control. (D) Analysis, by qRT-PCR of the relative expression of the 
afaABCDER genes from the AL51 1 AafaRr.KmFRT and AL51 1 strains grown to an OD 60 o of 0.6. Note: the afaE gene was naturally not expressed 
in the E. coli AL5 1 1 wild-type strain. 



afaD mRNA in the BW25113 + pILL1320 and 
BW25113 + pILL1323 strains (data not shown). The 
AfaR-dependent cleavage site at position A 66 within the 
afaD mRNA was confirmed, and an additional cleavage 
site was revealed at position G29 of the afaD mRNA 
(Figure 6F). 

RNA secondary structure prediction close to the RBS 
of afaD (Figure 2A) showed an AU-rich sequence next to 
a stem-loop helix and a G residue positioned two nucleo- 
tides upstream from a cleavage site mimicking the RNase 
E recognition sequence (47). Thus, cleavage in the 5' UTR 
and, probably, within the coding region of the afaD 
mRNA were dependent on RNase E. 



The expression of the afaR gene is controlled by 
temperature and a <r E promoter 

The AfaR gene may be transcribed by the alternative a E , a 
sigma factor associated with the response to extra- 
cytoplasmic stresses, such as increasing temperature. We 
therefore assessed dependence on the ct e promoter by 
introducing the pQFafaR::lacZ plasmid into E. coli 
wild-type and isogenic ArpoE mutant strains and assessing 
its activity by measuring P-galactosidase activity 
(Figure 7A and B). The afaR promoter had no transcrip- 
tional activity in the ArpoE strain, in either the exponen- 
tial or the stationary phase. A comparison of afaRr.lacZ 
fusion mRNA levels in the two strains confirmed these 
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Figure 5. Expression of the gene was required for the AfaR-de- 
pendent regulation of afaD mRNA decay. (A) Assessment of the 
relative expression of the genes of the afa-8 cluster in strain JW4130- 
1 (Ahfq) carrying pILL1320 versus BW25113 carrying pILL1320, 
determined by qRT-PCR. Asterisks indicate significant differences. 
(B) Assessment of the relative expression of the genes of the afa-8 
cluster in strain JW4130-1 (Ahfq) carrying pILL1322 versus JW4130- 
1 carrying pILL1320, determined by qRT-PCR. (C) Analysis by 
Western blotting and qRT-PCR of afaDr.gfp and AfaR gene expression 
in E. coli strain TOP10 Ahfq::FRT harbouring pZE2K-afaR or pZE2R- 
null plasmids together with pXG-0 (no target control) or 
pXGafaD::gfp target expression plasmids. The four isolates were 
cultured in LB medium at 37°C, to an OD 600 of 1.0. Expression 
levels of the gfp fusion gene were normalized to 1.0 for TOP 10 
Ahfqr.FRT carrying pZE2R-null and pXGafaDrgfp. AfaR expression 
was normalized to 1.0 for TOP10 Ahfqr.FRT carrying pZE2R-afaR 
and pXG-0. GroEL was used as a loading control. 



findings (Figure 7C). Thus, AfaR transcription is depend- 
ent on a E . 

Previous experiments had indicated that expression of 
the afa-8 gene cluster was temperature dependent and sug- 
gested that the a E -dependent AfaR sRNA was responsible 
for the thermoregulation of AfaD production (C. Le 
Bouguenec, personal communication). This hypothesis 
was tested by growing the BW25113 + pILL1320 strain 
at 20°C, 30°C, 37°C and 42°C, to an OD 600 of 0.6. The 
relative levels of the afaA and afaD mRNAs and AfaR 
were analysed by qRT-PCR (Figure 8 A), and AfaD pro- 
duction was assessed by Western blotting (Figure 8B). 
AfaR levels were higher at 37°C than at any other tem- 
perature. However, although afaD mRNA levels were 
high at 37°C, the level of AfaD production was lower, 
indicating that temperature controlled AfaR expression 
and the fine regulation of AfaD production. We tested 
this hypothesis by performing the same experiments with 
the BW251 13 +pILL1323 strain. We found that knockout 
of the afaR gene led to the constitutive production of afaD 
mRNA and AfaD (Figure 8C and D). 

The abundance of an sRNA at a given time point is a 
reflection of its rate of synthesis, mode of regulation and 
decay. We assessed the stability of the AfaR sRNAs by 
performing qRT-PCR analysis with total RNA isolated 
from the BW25113 + pZE2R-a/tfi? strain cultured at 
37°C and 42°C to an OD 600 of 0.4 before treatment with 
rifampicin to inhibit RNA synthesis. The half-life of AfaR 
was found to be 6.9 mini 21s at 37°C and 
5.04 min ± 17.7 s at 42°C. Given that the AfaR sRNA 
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Figure 6. RNase E dependence of the regulation of afaD expression by 
AfaR. (A) Schematic representation of the construction of the 
pILL1323 plasmid by allelic exchange from pILL1320. qRT-PCR 
analysis confirmed an absence of afaR expression from pILL1323 
(See Supplementary Figure SI). (B-D) Relative levels of the afaA and 
afaD mRNAs and AfaR sRNA, determined by qRT-PCR. in wild-type 
and RNase E ls thermosensitive strains. The values for each gene are 
expressed as a ratio of expression at 42° C (low level of RNase E) to 
that at 30° C (high level of RNase E) for the N3431 +pILL1320 strain 
(B, [RNase E ts ]), the N3433 + pILL1320 strain (C, [RNaseE + ]) and the 
N3431+pILL1322 strain (D, [RNase E ,s ]). (E) Circular RACE 
mapping of the afaC mRNA ends. Total RNA from E. coli 
BW25113 carrying pILL1320 (wild type) or pILL1323 (AafaR) was 
circularized by end-ligation with or without T4 RNA ligase (lanes + and 
— , respectively). The ligated 5' and 3' ends of the fragment were then 
amplified by RT-PCR. PCR products were analysed by electrophoresis 
in a 3% agarose gel. The band of interest was excised, cloned and 
sequenced. A ~310bp DNA fragment (gray arrow) was more 
abundant after ligation treatment (lane +) in the wild-type strain 
than in the AafaR strain, indicating an AfaR-dependent amplification 
of the 5' (cleavage at position -32) and 3' ends of the afaC mRNA. This 
suggests that AfaR is involved in the RNase-dependent cleavage of 
the afaABCD mRNA. (F) Location of RNase cleavage sites in the 
afaABCD mRNA. 
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was more abundant at 37°C than 42°C (Figure 8), the 
apparent slightly lower levels of sRNA at 42°C are not 
significantly different from those at 37°C (without mRNA 
target), providing no evidence of lower stability. The sta- 
bility of individual sRNAs varied widely, from highly 
stable to highly unstable, as reported for sRNAs in 
other bacteria (48,49). Taken together, these results 
indicate that AfaR transcription is controlled by 
temperature. 



DISCUSSION 

The invasion of a host by an ExPEC strain requires rapid 
adaptation to the new growth conditions associated with 
the transition from a free-living lifestyle to a host- 
associated state. This involves a reprogramming of the 
transcriptional activity of the bacterium. sRNAs 
function as gene regulators, relaying information from 
environmental sensors to response genes. We have previ- 
ously identified ExPEC-specific sRNA genes, including 
four sRNA genes within a pathogenicity island of the 
E. coli AL862 strain, close to gene clusters encoding 
known virulence factors (6). Here, we characterized 
SQ109, a candidate sRNA renamed AfaR, which is 
encoded by a gene located in the afa-8 gene cluster 
encoding a virulence-associated adhesive sheath. The a E - 
and temperature-dependent expression of AfaR is 
associated with the repression of AfaD-VIII invasin pro- 
duction, suggesting that AfaR is a stress-response regula- 
tor involved in the sensing of environmental stimuli 
relevant to host cell invasion. 

The synthesis of outer membrane proteins (OMPs) is 
controlled by a plethora of sRNAs, acting as antisense 
RNAs, in gram-negative bacteria (50-52). Some sRNAs 
are associated with OMPs involved in virulence. The 
InvR, RnaG and VrrA sRNAs have recently emerged as 
key regulators of the OmpD, IcsA and OmpA OMPs in 
pathogenic Shigella flexneri, S. typhimurium and 
V. cholerae, respectively (21,13,53). We previously 
showed that the FimR sRNA controls the production of 
type 1 fimbriae (an adhesion factor of the chaperone/usher 
family) involved in the colonization and invasion of the 
urinary tract (54). This sRNA acts by an antisense mech- 
anism, targeting the mRNA encoding the OMP FimD (6). 
In this study, we characterized AfaR, an sRNA that regu- 
lates afaD mRNA decay at the posttranscriptional level, 
thereby affecting AfaD-VIII invasin production. 
Recognition of the cellular pi -chain integrin by AfaD- 
VIII has been implicated in the internalization of beads 
coated with purified AfaD-VIII in cells derived from the 
uroepithelium, a target of a/a-5-carrying strains, suggest- 
ing that such strains are invasive (25). The AfaD-VIII 
invasin is not strictly an OMP, but our data suggest that 
it is associated with the bacterial outer membrane. The 
identification of AfaR thus enlarges the sRNA family of 
antisense RNAs regulating cell surface proteins involved 
in pathogenicity. 

During the biogenesis of adhesive structures of the 
chaperone/usher family at the cell surface, structural 
proteins are required in greater abundance than accessory 



proteins. Discoordinate gene expression is therefore essen- 
tial. Regulatory mechanisms involving the RNase-depend- 
ent and -independent processing of mRNA from E. coli 
gene clusters encoding P-, S-fimbriae, CFA/I and F1845 
adhesins have been reported (55-58). However, no sRNA 
has ever before been implicated in these processes. The 
identification of three sRNA gene candidates (fimR, 
afaR and prsR) overlapping adhesive structure-encoding 
gene clusters in a previous publication suggested a possible 
global involvement of sRNAs in the regulation of such 
virulence factors (6). In a putative sRNA-independent 
pathway of regulation, RNase E and RNase III have fre- 
quently been implicated in the processing of mRNAs for 
the biogenesis of adhesive structures (56,57). We therefore 
assessed the RNase dependence of the afa-8 gene cluster 
(Supplementary Figure S3A) and found that RNases E, 
III and LS but not RNases G and P, played a potential 
role in mRNA processing. Our findings demonstrate that 
the cleavage activity of RNase E in processing the 
afaABCD mRNA upstream from afaD is dependent on 
AfaR. Other RNase-dependent cleavage events remain 
to be fully characterized and did not seem to involve the 
AfaR sRNA directly. This suggests that two separate se- 
quential Afa-VIII biosynthesis regulation pathways may 
occur. The first involves cleavage of the afaABCD mRNA 
for regulation of the levels of each independent afa 
mRNA. In the second, the AfaR sRNA then specifically 
regulates afaD mRNA decay, probably switching AfaD- 
VIII invasin production on/off at the right time. A 
putative model of the regulation of the afa-8 gene cluster 
is presented in Supplementary Figure S3B. Furthermore, 
the afaD mRNA was found to be subject to several AfaR- 
dependent RNase E processing events. These cleavages 
induced rapid afaD mRNA decay, but low levels of trans- 
lation to generate AfaD protein nevertheless occurred. 
This suggests that AfaR is a finely tuned regulator that 
maintains AfaD protein synthesis in some conditions and 
facilitates invasion in certain circumstances. Thus, FimR 
and AfaR are the first sRNAs to be implicated in the 
regulation of OMPs exported by the chaperome/usher 
pathway. 

Pathogenic bacteria are highly adaptive organisms in 
which a tight control of virulence gene expression is essen- 
tial, and sRNAs are increasingly considered to be key 
elements of virulence control in E. coli. Our results 
indicate that production of the AfaD-VIII invasin is 
upregulated when AfaR is downregulated (Figure 8A 
and B), at temperatures above 37°C. The sRNA genes 
are, in general, regulated by various mechanisms unrelated 
to RNA stability, although sRNA regulation via tran- 
script stability has not been ruled out in most cases 
(48,49). Temperature does not seem to affect AfaR stabil- 
ity or AfaR decay in the absence of its mRNA target, 
indicating that the stability of the AfaR transcript may 
be crucial for the regulation of afaD mRNA decay at 
any temperature. First, temperature changes may inform 
the bacterium that it is within the host, triggering produc- 
tion of the AfaD-VIII invasin. Second, it might indicate 
that in vivo, variation of the host temperature should not 
affect the AfaR regulatory pathway. Lastly, the bacteria 
can produce more AfaD-VIII invasin at higher 
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Figure 7. a E is required for induction of the transcription of afaR. 
(A, B) P-Galactosidase assay on E. coli MC1061 (WT) and 
CAG22216 (ArpoE::Cm) strains carrying the pQF50 or 
pQ¥afaR::lacZ plasmid. The promoter of the afaR gene was dependent 
on a E sigma factor in the (A) exponential and (B) stationary (24 h of 
culture) phases. The absence of o E totally abolished expression of the 
afaRr.lacZ fusion gene. (C) Expression analysis of the afaRr.lacZ 
fusion mRNA relative to the 5S gene in stationary phase. AfaR expres- 
sion was abolished in the strain lacking a E . The asterisks indicate sig- 
nificant results. 



temperatures, presumably helping then to evade the host 
response to infection. Several putative temperature 
sensors might detect and relay information to the AfaR 
sRNA. We tested the hypothesis that the a E factor plays a 
crucial role in this process. Indeed, a E is a sigma factor 
associated with the extracytoplasmic/heat stress responses 
(59), as observed for expression of the P-fimbria gene 
cluster, which is controlled by temperature in a a -de- 
pendent manner (60). More generally, the alternative a E 
factor upregulates the expression of many other sRNA 
genes, including MicA and RybB in Salmonella species 
(61) and, probably, FimR in ExPEC strains (6). Our 
data provide new evidence that the a E regulon may play 
a more important role in the virulence of gram-negative 
bacteria than previously suspected. 

Expression of the afa-8 gene cluster may also be 
controlled by the AfaF and AfaA proteins, which are 
homologues of the transcriptional regulatory proteins of 
the PapI/PapB family. We identified, in silico, six putative 
AfaA-binding sites (5'-CTATATTTT-3' consensus 
sequence) in the afaA and afaF intergenic region and 
three other binding sites overlapping the promoter 
region of the afaR gene (Figure 4A). These findings 
strongly suggest that afaR transcription may also be de- 
pendent on these regulatory proteins. We also found that 
AfaR expression was strongly downregulated at 42°C, 
suggesting that AfaR transcription may be dependent on 
another, as yet undetermined mechanism. Data from 
several studies have shown the CpxAR two-component 
systems control the production of curli and P-fimbriae in 
uropathogenic E. coli, of type IV bundle-forming pili in 
enteropathogenic E. coli, and downregulate the expression 
of CT E -dependent genes, rendering them of interest for 
future studies (60). 

Comparative genomics of the afaDjafaE intercistronic 
region of the afa-8 gene cluster showed it to be very 
specific in size and composition, and no sequence related 
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Figure 8. Control, by temperature, of the production of the AfaD 
invasin by AfaR regulation. (A) Assessment of the relative levels of 
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grown to an OD 600 of 0.6 at 30°C, 37°C and 42°C versus 20°C, as 
determined by qRT-PCR. The expression levels of all genes at 20°C 
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tion by E. coli BW25115 carrying pILL1320 in the same growth con- 
ditions as (A). (C, D) as in (A, B), but with the BW25113 strain 
carrying pILL1323 (AafaR: : KmFRT). GroEL was used as a loading 
control. 



to the afaR gene was detected in the other gene clusters of 
the afa family. Thus, AfaR seems to be specific to the afa-8 
gene cluster. We also found that the 3' end of the afaR 
gene displayed sequence identity to an IS2 insertion 
sequence, at both the DNA and protein levels. There is 
therefore evidence to suggest that afa-8 evolved from its 
ancestor at least partly through IS2-mediated disruption 
of the afaDjafaE intercistronic region. Pseudogenization 
of the IS2 transposase gene may then have led to the 
creation of the afaR gene. This reorganization of the 
intercistronic sequence also seems to have had conse- 
quences for the transcription units specific to the afa-8 
gene cluster. Indeed, only in afa-8 was it was possible 
to identify an afaE promoter and a rho-independent ter- 
minator for the afaABCD mRNA. These promoter 
and terminator sequences may have been acquired 
from the IS2 sequence, as previously reported follow- 
ing the insertion of IS2 into the ampC promoter (62) 
or of IS3 insertion into the 5' end of argE (63). IS 
elements have been shown to play a role in the dispersal 
of SprA sRNA genes between strains of the genus 
Staphylococcus (9). Finally, our study suggests that 
some of the sRNA genes in E. coli may be derived 
from IS elements. 
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